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Abstract - We consider the spectrum of the evolution operator for bound 
chaotic systems by evaluating its trace. This trace is known to approach 
unity as t ^ oo for bound systems. It is written as the Fourier transform of 
the logaritmic derivative of a zeta function whose zeros are identified with 
the eigenvalues of the operator. Zeta functions are expected to be entire 
only for very specific systems, like Axiom-A systems. For bound chaotic 
systems complications arise due to e.g. intermittency and non completeness 
of the symbolic dynamics. For bound intermittent system an approximation 
of the zeta function is derived. It is argued that bound systems with long 
time tails have branch cuts in the zeta function and traces approaching unity 
as a powerlaw. Another feature of bound chaotic systems, relevant for the 
asymptotics of the trace, is that the dominant time scale can be much longer 
than the period of the shortest periodic orbit. This the case for e.g. the hy- 
perbola billiard. Isolated zeros of the zeta function for the hyperbola billiard 
are evaluated by means of a cycle expansion. Crucial for the success of this 
approach is the identification of a sequence of periodic orbit responsible for 
a logaritmic branch cut in the zeta function. Semiclassical implications are 
discussed at the end. 
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1 Introduction 



Time correlations provide an efficient tool for the analysis of dynamical systems. For, 
some, chaotic systems correlations decay exponentially, reflecting loss of memory during 
evolution. Often the correlation function is found to be modulated by some typical 
frequence(s) related to the periodic orbits of the system. Its Fourier transform will 
thus have poles {resonances) in the complex frequence (/c) plane. The imaginary parts 
correspond to the decay rates and the real parts to the oscillation frequences. Integrable 
systems have resonances on the real axis, reflecting quasiperiodicity of the motion. 

A rigorous theory for resonances has so far been developed only for Axiom-A systems 
|]l|, 1^, 1^. The correlation function is then meromorphic (i.e. it has nothing worse than 
poles) in a strip around the real A;— axis. For mixing Axiom-A maps one can establish 
exponential bounds on the decay of correlations, i.e. there is a strip (gap) above the real 
axis without resonances. But even for Axiom-A flows the resonances may lie arbitrarily 
close to the real axis leading to non-exponential decay of correlations |5|, ^. It is a 
general experience that the positions of the resonances are independent of the particular 
observables whereas the residues are not. This is proven only for Axiom-A systems but 
numerical work suggest that it holds more generally. This is very intriguing since one can 
thus characterize a system by its (classical) resonance spectrum, in nice analogy with its 
quantum spectrum. 

However, Axiom A systems are very idealistic systems from a physical point of view, 
although they may be realized by e.g. some open scattering systems. The purpose of the 
present paper is to study resonance spectra in bound Hamiltonian systems, exhibiting 
chaos. Then we know that there is a trivial pole at /c = reflecting the boundedness 
of the system. All other resonances lie in the upper halfplane. This is essentially all we 
know for certain. 

Crucial for the succesful theory of Axiom-A system is the existence of a finite markov 
partition. One can thus perform a well converging cycle expansion of a zeta function 
in order to locate the resonances Such finite grammar does generally not exist for 
bound systems. This leaves us without any a priori knowledge of the analytical structure 
of the zeta function we are going to study. 

We will restrict ourselves to ergodic billiards in the Euclidean plane. Seemingly all 
systems of this type have one thing in common; they are intermittent. Typically, the 
motion will intermittently alternate between a chaotic region of phase space and a regular 
one. This might imply power law decay of correlations. In a paper by Baladi, Eckmann 
and Ruelle Q (hereafter referred to as the BER approximation) it is demonstrated how 
the approximate positions of the resonances for intermittent systems may be calculated. 
This result is directly applicable to a wide class of bound chaotic systems and we will 
make use of this observation in the following. 

The second main theme of the paper is a periodic orbit sum rule discovered in It 
is thus shown that a certain weighted sum over all periodic orbits up to period t tends 
to unity as t ^ oo. This is an asymptotic result and it is of course highly interesting to 
know how this asymptotic limit is approached. The sum rule is a crucial ingredient when 
computing correlations in quantum spectra. 

The paper is outlined as follows. We begin, in section 2.1, by reviewing the derivation 
of the classical zeta function whose zeros are associated with the resonances. Then, in 
section 2.2, we review how, according to the BER approximation, the positions of the 
resonances are related to a simple probability distribution associated with the system. 
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In section 3 we relate, by heuristic arguments, different features of this probabihty dis- 
tribution to the resonance spectrum. In section 4 we apply cycle expansions of the zeta 
function to a specific dynamical system - the hyperbola billiard. In section 5 we discuss 
the relevance of our result for correlations in quantum spectra. 

2 Theory 

There are mainly three methods for the determination of resonance spectra. The pedes- 
trian way is to numerically simulate the system, compute the correlation function for 
some physical observable(s), make a Fourier transform and then locate its poles. This 
is of course very tedious and the success moderate and the theoretical understanding 
gained minor. We will thus focuse on the remaining two methods. 

2.1 The zeta function 

Periodic orbit expansions (or cycle expansions - we will use the words cycles, periodic 
orbits or simply orbits synonymously in the following) of zeta functions have shown 
successful for calculating various averages of chaotic sets. Most of the labour has so far 
been devoted to maps, that is, systems with discrete time. A zeta function for systems 
with continous time has been derived in [^]. We briefly review the derivation. The 
resonances /c„ may formally be associated with the eigenvalues exp{iknt) of the evolution 
operator £*, acting on a function as 

£'^{x) = J 5(x-f\y))^{y)dy . (1) 

The phase space point x is taken by the flow to /*(x) during time t. Our main concern 
in this paper is to compute the trace of this operator, that is, the sum of its eigenvalues. 
The trace may be written as a sum over the isolated periodic orbits in the system 

where n is the number of repetitions of primitive orbit having period Tp and Mp is 
the Jacobian of the Poincare map. This is exactly the sum occuring in ref. |^] which is 
shown to approach unity as t ^ oo, provided the system is bound (and no other system 
will be discussed from now on). 

The function tr{t), (we will refer to this object as simply the trace in the following) 
may be written as the Fourier transform of the logaritmic derivative of a zeta function: 



We restrict ourselves to systems with two degree's of freedom for which all periodic orbits 
are isolated and unstable. The zeta function then reads 



where Ap is the expanding eigenvalue of Mp. Each zero ka of Z(k) induce a pole in the 
integrand and is identified with an eigenvalue of the evolution operator. The leading zero 
is the escape rate, which is zero for a bound system; ko = 0. 

The corresponding zeta function for systems with a hyperbolic analytic Poincare map 
has been shown to be entire provided the existence of a finite Markov partition If 
the zeta function is entire we can write 



tr{t) = l + J2 (^) 



In eq. (|5|) we have neglected the integral along the upper semicircle. It may be argued 
that this can only yield a deltafunction 6{t) and derivatives of the deltafunction 

If there is a gap between the real axis and the nontrivial zeros then the trace approach 
unity exponentially fast. 

A an example consider the following simplified Axiom-A system (it may be realized 
by a three disk scatterer with the disks infinitely separated) . We take the system to have 
a binary complete symbolic dynamics. An orbit p has eigenvalue Ap = Aq'' and period 
Tp = ripTo where Up is the length of the symbol code p. Expanding the product over 



cycles gives, c.f. ref. |10| 



^(^)=n(i-2^) • (6) 

m=0 \ ^^0 / 

The product is convergent so we can just read of the zeros 

kn,m = ^ {27m + i[{m + I) In Ao -In 2]) n = 0,±l±2... . (7) 

The spacing between zeros in the real direction is directly related to the single timescale 
in the system - the period of the shortest periodic orbit. The situation for bound systems 
is generally much more complicated, as we will see later; they may exhibit an infinite set 
of time scales. 

For a more general system with complete binary grammar the spectrum above will 
give the gross structure, whereas longer periodic orbit provide corrections |1C, How- 



-ikTo ^ ™+l 



ever, even for small deviations from the simple example above there will be quite intricate 
cancellations between poles and zeros [Q, |l2| of the individual m-factors. For bound 
systems there is no finite Markov partition and no hyperbolicity. There is then no reason 
to expect entireness of the zeta function and we will concentrate on revealing the main 
singularities. 



2.2 The BER approximation - a probabalistic approach 

In ref. a method has been given for the determination of the approximate resonance 
spectrum for intermittent systems. In such a system there are two, more or less, distinct 
phases; one regular and one irregular (chaotic). Call the consecutive instants when the 
system enters the regular phase {ti} and consider the intervals /j = Provided 
the chaotic phase is chaotic enough, the motions in different intervals are nearly mutually 
independent. In particular the lengths of these intervals Aj = ti — tj-i are mutually inde- 
pendent and A may be considered as a stochastic variable with probability distribution 
p(A). Consider now the Fourier transform 

p(k) = / e-^^^p(A)dA . (8) 
Jo 
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Under the assumption above, the roots of 

1 - p{k) = (9) 

will now provide the approximate resonances. We will modify the calculation of ref. 
so as to apply directly to the trace of the evolution operator. This will enable us to find 
a very simple approximation of Z{k) in terms of p(A). 

First we turn the phase space integral in eq. (§) into a time average. 

tr{t) = V{5{x{to)-xito + t)))to . (10) 

where V is the available phase space volume. 

Let now trmit) be the trace conditioned by 'if to ^ In then t + to € In+m ■ The 
trace is now expressed as the sum tr{t) = X]m=o ^'^m (^) • We now define two probability 
distributions: 

• u) is the probability that the system is at phase space point x at time to and 
that there is time u left of the current interval. 

• f2{y, v) is the probability that the system is at point y at time to + t and that time 
V has ellapsed since the current interval was entered. 

We can now express the average as 
tr^{t)^V <A> [ dudvp^'^-^^t-u-v) [ dxdy5{x-y)fi{x,u)f2iy,v) , (11) 



where p*"(t) is the n-fold convolution of p(t), and < A > the mean length of the intervals. 
We now express the convolution by means of the Fourier transforms 

trmit) « ^ ^ / dk e'^^p{kr~^ [ du dv dx dy 5{x - y)e*'="/i(x, n)e*^V2(y, t^) = 



27r 



y < A > 



27r 

We can now sum over m and arrive at 



dk e'''^p{k)"^-^ I dxfi{x,k)f2ix,k) .(12) 



tr(t) « tro(t) + ^^^^^ I dk e*'*Y37(fc) ' / dxhix,k)hix,k) . (13) 



dxf\{x,k)f2{x,k) . (14) 



Comparing with eq. (^) we indentify 

Z'jk) _ V <A>i 
Z{k) ~ l-p{k) 

Moreover, we can identify: 

Z{k) ^ Z{k) = l-p{k) (15) 

provided that 

V <A> J dxfi{x,k)f2{x,k) = j Ae-'^^p{A)dA . (16) 

To verify that this is indeed the case we make use of the fact the remaining time of the 
current interval n is a function of the phase space point u = gi{x) so that 

fi{x,u) = <p{x)6{u - gi{x)) (17) 



5 



where (/'(x) is the phase space density. Similarly we have 

f2{x,v)=(l){x)6{v-g2{x)) . (18) 

We can now express the probability function p(A) weighted with the length of the interval 
A in terms of gi and g2 

A • p{A) =<A> J 0(x)5(A - gi{x) - g2{x))dx (19) 

Also using the fact that the phase space density is uniform: (l){x)-V = 1 eq. ( p^ ) is easily 
verified. 

In this section we have been able to give the zeta function a probabalistic interpre- 
tation and have constructed an approximate zeta function Z{k) whose zeros should be 
interpreted as the approximate resonances. Both functions, Z{k) and Z{k) have the 
trivial zero feo = for unitarity reasons. 

3 Phenomenology of the BER approximation 



The close relationship between Z{k) and Z{k) as defined by eq. (|T5|) motivates a study 
of how different features of the probability function j?(A) affects the spectrum and trace 
of the evolution operator. Throughout this section we will ignore the difference between 
Z{h) and Z{k), that is we will temporarly assume that Z{k) = Z(k). The reader should 
keep in mind that the spectra and traces obtained here is of course only approximations 
of those of real dynamical systems and the errors are not under control. 

We begin by studying a model (model A) with a sharp step and and exponential 
decay 

P.(A) = { l^..,r, i ^ . (^») 

where is a normalization factor. We take the Fourier transform and get 

1 + ikT2 

One can obtain the following asymptotic expression for the zeros of this function 

fc„p^^(±27r(n-i) + i-ln(f 27r(n-i))) n = 1,2,3,... , (22) 

which agrees very well for large n and/or large c.f. fig. la. We see that the relevant 
time scale is Ti, whereas T2 is of minor significance. Z^ik) has a simple pole at k = i/T2 
reflecting the behaviour of pyi(A) at infinity. This is all we need to calculate the trace 
by means of eq. (|^. We get 

tr{t) = ^,eros{t) + '^tail{t) , (23) 

where 

^™(t) = l + E^^'"* « l + (-^)-*/^^ 2cos(2vr(n--)-)-(2vr(n--))-*/^^ , (24) 

n n=l 1 
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and ^ 

^taii(i) = -e"^ • (25) 
The sum ( p^ ) is absolutely convergent when ^ > 1 and tends exponentially to unity 
as ^ ^ oo. To study its behaviour for small jr we have to the smear the function term 
by term, that is convolute it with a gaussian having width a. For a general spectrum 
kn = ior, + i • 6„ we find for the smeared sum 



^2 



'^zeros,a{t) = 2 ^ cxp |^ — —{al - b^) - bntj cos{ant - a'^anbn) ■ (26) 

For all relevant spectra this sum is absolutely convergent for any smearing width a. 

The trace when ^ = 0.5, is plotted in fig. lb. The sharp step at ^ = 1 reflects 
the sharp step in p^(A). The are indeed steps whenever jr is an integer, but due to the 
exponential damping they are less pronounced. 

Let us now consider a case with power law decay of p(A) (model B). 

p^^^^ = [ Ni^r A > . 

Again we compute the Fourier transform and seek the zeros of 

ZB{k) = {m-l)-Em{ikTi)-l , (28) 
where Em{z) is the exponential integral. It turns out that all zeros a positioned where 



we can use the asymptotic expansion of Em{z), see [13]. The position of the zeros are 
well approximated by 

fc„«^(^±27r(n-i)+i.ln(H!g^)^ n=l,2,3,... , (29) 

which is equivalent with model A when m — l = ^ , that is when the height of the step in 
p(A) is the same. High zeros (large n) are thus only sensitive to the immediate vicinity of 
the step. In fig. la the exact zeros of model A and B are compared with this asymptotic 
result. Formally, if we increase m (or ^ in model A) ^^(A) will tend to a delta peak, 
and the resonances approach the real axis, resembling an integrable system. Of course, 
the assumption in the BER approximation is then no longer fulfilled. 

However there is one major difference to model A. The zeta function is no longer 
meromorphic since Em{ikTi) has a branch cut along the positive imaginary A;— axis. 
This is a consequence of the power law tail of pb(A). 

We therefore have to modify the contour used in section 2.1 

where ^branch is given in fig. 2. We now get 

tr{t) = ^,eros{t) + ^tail{t) , (31) 

with 

^^tauit) = r -e-^^Im{ ^'l'f^'h dC , (32) 



7 



and ^ 

zerosii) before, e is a small positive number. The resulting trace for model B is 
shown in fig. lb. The discrete ^zerosit) is the same as in the previous model (A) since 
we are using the asymptotic expressions (|2^) and (2£). It is straight forward to show 



that ^taii{t) —1; t ^ and ^taiiit) —>■ t —>■ oo. The asymptotic behaviour of 

the trace for general p{A) with power law tail is given at the end of this section. 

Let us modify the first model (A) so as to have a discontinuity instead of a C". 

P^^^^ = [ iVexp(I^) A>ri . (33) 

There is again a pole at k = i/T2. The zeros are now approximated by 

fcn ~ fT- (±27rn + 2iln(g27rn)) n = 1,2,3,... . (34) 

Apart from shift a in the real direction, as compared to model A (and B), the important 
effect is the factor 2 in the imaginary part - the zeros lie further away from the real axis. 
The corresponding trace is given in fig. 3. 

Finally let us demonstrate how small irregularities in p(A) will have a drastic effect 
on the spectrum. We will consider 

f A < T 

Pd{A) = 1 iV(A-r)exp(-|) r<A<2T (35) 

[ iV(A-r)exp(-|)(l + e) 2r < A . 

Z£)[k) has a double pole at A; = i/T. The spectrum of zeros for different perturbations, e, 
is given in fig. 4b. In the low part the unperturbed spectrum is only sligthly affected. The 
perturbation corresponds to a nonleading generations of resonances. But eventually the 
two generations interfere, leading to a modulation in the spectrum. For high t the higher 
time scale (2T) takes over and dominates. This result leaves us with considerable doubt 
on the possibility on determining the resonance spectrum from a numerically obtained 
probability distribution p(A). If there is a single dominating time scale it seems perfectly 
feasible [0]. But there might hide further time scales, with considerable potential impact 
on the spectrum. The traces for e = and e = 0.01 are given in fig. 4c. For the 
unperturbed case, the trace is essentially zero when < A < T. This is a remarkable 
conspiracy between ^ zeros and ^taii = ~2 • e~*/^ to cancel in this way. But from a 
physical point of view it seems sensible, c.f. eq. (§). There is, in any system, a shortest 
periodic orbit, with period Tmin, such that tr{t) = when t < Tmin- 

When we add a perturbation, e, the trace changes in two respects. First, pd{A) has 
a small step at t = 2T, which is not visible in the figure. When < A < T the trace is 
no longer positive definite. This is not to allarming since we have neglected trQ{t) in eq. 
(p!3[). Moreover, the spectral sum is divergent in this region so it is questionable what 
confidence we should have in the results there. For larger t we see that a small change 
of p{A) corresponds a small change in the trace. 

Let us summarize our findings so far, remembering that our evidence is only heuris- 
tical and that we have assumed the exact equality Z{k) = Z{k). The zeros of the zeta 
function depends almost entirely on the breakpoints (i.e. different kinds of discontinu- 
ities) in the probability distribution p(A). Each such breakpoint introduces a time scale. 
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When several such timescales are present the largest will dominate sufficiently high up in 
the spectrum, in the sense that the spacing between succesive zeros in the real direction 
is directly related to this timescale. 

The zeta function is entire only if p( A) decays faster than exponentially. If the decay 
is exponential the zeta function will have a pole on the imaginary axis. In this case the 
trace will approach unity exponentially fast. If p(A) decays as a powerlaw there will be a 
logaritmic branch cut along the positive imaginary axis implying powerlaw convergence 
of the trace. It is straightforward to show that if p{A) ~ l/A*", the tail contribution, 
"^taiiit), to the trace, will behave asymptotically as 



AqXAq 

(m - 1) < A > t*^ 



^^taiiit) - , -,1^ 00 (36) 



with m >3, where <> denotes mean value. Aq is any point in the tail of p(A). 



4 Cycle expansions - studies of a dynamical model 

In this section we will study periodic orbit expansions of the zeta function @) for a 
specific dynamical system, namely the hyperbola billiard, i.e., a point particle elastically 
bouncing off the walls given by the equation 

\xy\ = l , (37) 



An effective algoritm exist for the computation of the cycles in this system For 
a definition and discussion of the symbolic coding of cycles used in this paper, please 



consult ref. |15]. We use units such that the particle mass m = 1 and we fix the 
energy to E = ^, the length and time scales thus coincide. The zeta function may be 
factorized into the irreducible representation of the symmetry group which is 
for the hyperbola billiard. We will restrict our attention to the symmetric representation 



Ai containing the trivial zero /cq = 0, |16|. This is equivalent to restricting the system 
to the fundamental domain |15, 16 1. The periodic orbits considered in the following are 
periodic in this fundamental domain. They have to be retraced 1,2 or 4 times to close 
onto themselves in the full system. Although we will frequently draw pictures of the 
periodic orbits in the full system since they are easier to visualize that way. 
First we expand the m-product of @ 

p m=0 \ I PI p / p np=0 I PI 

Reccurence relations for the expansion coefficients 6„ are derived in the appendix. If we 
now expand the p-product, we obtain a Dirichlet series consisting of all distinct linear 
combinations (pseudo orbits) = [rip],: 

Z{E) = J2CNe-'^" , (39) 
N 

where we have defined the quantities 
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TN = Y,n,Tp . (41) 

p 

There are no obvious choice of how to order the terms, but natural choices are to 
order them according to i) length of the symbol code, ii) increasing period T/v or, iii) 
decreasing amplitudes Cn- Our default option will be the last one. Switching between 
the three options will mean a complete reorganisation of the sequence the pseudo orbits. 
Our default choice means that we can make no simple statements about the conditional 
convergence of the expansion. This will however not bother us to much. Typically, 
nontrivial zeros of zeta functions are positioned where the products, or series, are not even 
conditionally convergent, c.f. the Riemann zeta function, although in some lucky cases, 
due to alternating signs of the Cat's, they will lie in the strip of conditional convergence 

The idea is now to compute the partial sums (the first Ncut terms) and determine 
the zeros of the truncated expansion and study their possible convergence when Ncut is 
increased. To this end we have computed all prime orbits with eigenvalue Ap < 2700, 
and all orbits with T < 14, both samples containing about 10^ orbits. To have something 
to crosscheck our result with, we calculate the trace by means of eq. (||). Of course, we 
have to use smearing so we replace the delta function in eq. (§) with a gaussian having 
width a = 0.5 The result is given in fig. 10. Due to the exponential proliferation of cycles 
this method soon becomes completely unrealistic as t is increased. 

Our first attempt is to simply include all orbits up to the cutoff. The zeros thus found 
will show no tendency of convergence (the zeros drift monotonously towards the origin as 
Ncut is increased), and the trace, computed by means of the sum over zeros (again using 
smearing, c.f. eq ([2^)), will not show the slightest resemblance with the exact trace. 
Such failure of convergence are usually due to some singularity in the zeta function. In 
order to identify this we will, in some detail, study the structure of the set of cycles. 
Inspired by the BER approximation we divide (somewhat ambigously) the billiard into 
one chaotic part in the center and one (almost) regular out in the arms. We can thus 
divide the cycles into three subsets. 

1. Cycles with bounces exclusively in the chaotic region. The eigenvalues A are ex- 
ponentially bounded with length. 

2. Cycles with bounces in both regions. 

3. Cycles with bounces exclusively in the regular arms. 

In the following we will focuse our attention on the types of cycles which dominate 
the expansion, with ordering chosen as above. These are the pseudo orbits with large Cn 
which means that the important cycles come from set 2 and 3 above - cycles spending 
most of the time in the arms. 

It is also convenient to distinguish between simple cycles, that is cycles making only 
one (or no) excursion into an arm, and composed, that is cycles making more than one 
such excursion. The simple cycles occur in sequences of the form pO^ where the row of 
O's corresponds to a number (~ j) of oscillating bounces in an arm and p corresponds to 
the behaviour in the central region. 

Some examples of simple cycles are given in the table and fig. 5. Cycles belonging 
to the sequences 20-^ and llO-' make one excursion into an arm and are then injected 
into the opposite arm without any bouncing in the central chaotic region. These two 
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sequences are in fact unbounded. All other simple sequences are, we believe, finite. A 
cycle 210-' makes a number of oscillations in one arm, sneaks into a neighbouring arm 
without bouncing in the middle. This is obviously a very short sequence. A cycle in 
the sequence 22(P makes a number of oscillations in one arm, makes one scattering in 
the central region and is then injected back into the arm it came from. A cycle 221(P 
makes a number of oscillations in one arm, bounces once in the central region and is then 
injected into a neighbouring arm. 

It is obvious from fig. 5 why finite families pO^ are terminated at some j = jmax- 
There is a section of the orbit in the central region coming closer and closer to a branch 
of the boundary as j is increased. When j = jmax + 1 the orbit would need to go through 
this branch which it of course is not allowed to. For this reason the finite sequences are 
naturally organised into families where one member go directly to the opposite branch 
when another member hit a branch on the way, see fig. 5a . Thus 210-' may be considered 
as a member of (210^10^+3). Other important familes are (220^2110^1210^1010^+^) 
and (2210^2010^■+^120J■+^ 1110^+2)^ Ah members of a family terminate at the same 
J = jmax, cf- refs 1 15, O]. The family member with the fewest bounces in the central 



region has the largest weight and the rest are of minor importance. 

The class consisting of infinite simple sequences, will be named Cmf in the following 
and contains only 20-' and llO-'.The class of simple finite sequences is named The 
orbits in these classes are taken from sets 2) and 3) above respectively. The orbits in 
subset 1) may be absorbed into e.g. class C fin- 
Composed cycles doing two excursions into the potential arms may be written on 
the form pi0^^p20^'^. Picking the simple cycles (or rather symbol strings) piO^^ and p20^^ 
from Cinf and/or Cjin we can construct three classes of such composed cycles: Cinf*Cinf, 
Cinf *Cfin and C/j„ *Cfin- Generalisation to higher orders are straightforward. Members 
of the 2-composed sequences 2(P^20^\ 20^1210^2 and 210^1210^2 (f^om the three classes 
respectively) are displayed in fig. 6. The latter two are again finite. The allowed cycles 
may be represented in the ji , j2 plane as in fig. 7. 

The region of allowed cycles is quite different for the class Cj„,/*Cj„,/. Take for example 
the 2(y^2(P'^ sequence. In the adiabatic approximation one can show that, in the limit 
ji ^ 1 and j2 3> 1, the allowed cycles lie in the infinite region 

-<j2<a-ji , (42) 
a 

where a 1158. 

Of course, the division between simple and composed cycles is again somewhat am- 
bigous; e.g. should 21210-^ be considered as simple or to belong to the composed sequence 
210{ 210^2 > c.f. fig. 5b and 6a? 

Let us now discuss the contributions of the different sequences to the expansion of 
the zeta function. 



We begin with the class Cinf- For large j it was shown in [19| that the periods and 
the stability eigenvalues for the 20-' and llO-' sequences are: 



Aj-c{j + l) j = l,2,...,oo , ^ ' 

where c = 117/4 for 20-'. The corresponding factor for llO-' is much bigger. We suspect 
these infinitite sequences to be the main cause of our divergence problem. Before following 
this track we discuss the finite sequences. 
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A finite simple sequence pO^ will have invariants depending on j similarly to the 
infinite ones above, but with j terminating at some jmax, see the table. The slope of Aj 
versus j will depend on the behaviour in the central region and seems to be exponentially 
bounded with the number of bounces there. Each finite sequence naturally introduces a 
time scale given by the period of the last orbit, some of these timescales are named in the 
table. The trace, fig 10, will obviously show a sudden decrease, cf. eq. (§) at any such 
time scale. In fig. 10 they are rather recognized as smooth peaks due to the smearing. 
These (infinite) set of timescales are essential for the problem and will be reflected by 
the zeros of the zeta function. 

The crossterms between simple orbits are naturally discussed on the same footing as 
the composed sequences. There is a potential possibility of shadowing, that is, the term 
due to orbit pi0^^p2(P'^ is compensated by the pseudo orbit consisting of pi(P^ and P20^'^. 
This can however only be partly realised due to the finiteness (pruning) of the sequences. 
The double hatched area in fig 7a shows where there exist compensating crossterms. 
However these terms are sparse in the expansion, and may in some sense be considered 
as corrections, cf. the curvature corrections in ref. |10|. The orbits along the boundary 



of the allowed region, cf. fig. 7, also introduce timescales. 

We suspected the class Cmf to be the main cause of our problem. We will now 
attempt a factorization 

Z{k) = Zi{k) ■ Z2{k) , (44) 

where, in the first factor, we only include cycles p G Cinf and in the second we include 
the rest p G Ctot — Cinf ■ The first factor is approximately 

ziik)^i--y2^— ■ (45) 

i=i 

Crossterms are sparse in the series, so we have omitted them. It is straight forward to 
see that the function defined by this sum have a logaritmic branch cut along the positive 
imaginary axis. 

Our hope is now that the eigenvalues are given by the zeros of Z2{k), which is ex- 
panded as before. We now seem to obtain a well defined spectrum where when Ncut > 180 
see fig. 8! However, we do not find much improvement when Ncut is increased further. 
It is thus obvious that the zeta function needs further regularization. The convergence 
of the trivial zero ko is plotted in fig. 9. 

We obviously need some support for our claim that the obtained spectrum is indeed 
close to the exact one. To that end we compute the trace by means of the sum over zeros 
^zerosit) = 1 + Z^n^*^"* ^^^^ compare with the exact result due to eq. @). Again we 
use gaussian smearing, which is unavoidable since we only have a finite sample of zeros 
of Z{k). The result is plotted in fig. 10. We see that the curves closely resembles each 
other on the small scale structure. 

Let us first study the structure of the spectrum, cf. fig. 8. We recognize several 
features from the spectra in section 3. The imaginary part of the zeros increase slower 
than logaritmically with the real parts (cf. model D in section 3). We also recognize a 
modulation of the spectrum. These features indicate that several time scales are involved. 
The real part of the zeros are more or less equally spaced. The mean spacing correspond 
to a time scale T ~ 27.8 which is much longer than the shortest periodic orbit of the 
system, Tmin = 2. It is indeed close to the timescale Tc ~ 27.61 associated with the 
sequence 220-' which dominates the expansion of Z2 for the first 180 terms. Clearly, 
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using ordering option ii) we would need all cycles with Tp < Tc in order to resolve the 
spectrum. This would mean ~ 10^ cycles. With the ordering option chosen we have 
managed to obtain a good approximation to the spectrum using a thousand times fewer 
orbits! 

In fig. 10 we see that the exact trace and the sum ^ zeros differ by a slowly varying 
function. The immediate suggestion to explain this discrepancy is that there is another 
contribution to the trace from a pole or a branch cut, c.f. ^taii- The latter possibility is 
suggested by the function Zi{k) having a branch cut. However we do not yet claim that 
Z2{k), as defined above, is entire. We cannot expect the branch cut to be factored out 
that easy. This means that we cannot compute this extra contribution only from Zi. 

The problems encountered so far could be somewhat enlightened by explicit knowl- 
edge of the probability distribution p{A) for the billiard since we expect there to be a 
close relationship between its Fourier transform and the zeta function. This is easily 
obtained by simulating the system. It is thus essential that we take advantage of the 
adiabatic invariance in the arms, c.f. ||l9|. There is no sharp borderline between the 
chaotic and regular part of the billiard. We make an arbitrary division at some x = Xdiv 
(remember the system is defined on the fundamental domain y > 0,y < x). The result 
is given in fig. 11. The histogram is based on 10® visits in the arm and the borderline is 
set to Xdiv = 4. 

We recognize the timescales as introduced above and defined in the table, as under- 
lying the structure of p(A). 

For small A < 15 the result is sensitive to the exact location of the borderline. But 
the rest of p(A) is rather insensitive to it. 

The curve exhibit a peak around A « T;,. Among the important composed sequences 
210*1210*2 and 20^1210^2 (^^th ji ^ 0) all cycles have T < T^, and indeed their periods 
are distributed around the peak of p{A). 

When Tc < A < Td, p{A) is decreasing exponentially to a high degree of accuracy. In 
the region of exponential decay there seems to be a gap, with no pronounced timescales 
(as defined by means of finite sequences). 

When A > Td, the curve slowly begins to deviate from the exponential. Relevant 
sequences for this tail are e.g. 20^120^^ _ f^j-^g sequence 2020-^ terminates at j ~ 990, 
T « 116. 

It would be highly desirable to know the asymptotic form of p(A), if possible we 
could be able to estimate the branch cut contribution ^taiiit)- However, the crossover 
from exponential decay to something else (powerlaw?) warns us that p{A) is not easily 
extrapolated. There may hide more surprises higher up. Estimating the tail behaviour 
by analytical means would require a detailed understanding of how the almost integrable 
arms are coupled through the chaotic central region. The problem lies in the very long 
laminar phases (large A) which have segments of the trajectory coming from one arm and 
injected directly into the opposite without bouncing in the chaotic part (in the simulation 
we require at least one central bounce to consider the laminar phase as terminated). From 
a cycling point of view the classes Cmf, Cm/ * Cinf etc. are central for the description 
of these trajectory. Excluding Cinf as a first step in the regularization, as we did, seems 
quite reasonable. However, to proceed the regularization scheme further will require 
entirely new ideas. Such a procedure would probably reveal more singularities. 

In the calculation of the spectra above we concluded that the timescale Tc is in some 
sense dominant. An explanation of this is offered by the shape of p(A). Tc is the last 
timescale (breakpoint) before the exponential part of p{A). We conlcuded in section 3 
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that such a behaviour of p(A) would lead to a mean spacing corresponding to Tc, at least 
in the lower part of the spectrum. We stress again that it is not in practice possible to 
calculate the spectrum from p(A) obtained by numerical simulation since it is the fine 
structure of p{A) that determines the spectrum and this may not be resolved by any 
reasonable statistics. 

Allthough this section ended in some frustration we summarize our main achieve- 
ments: We have computed several isolated zeros of the zeta function for the hyperbola 
billard. We have also presented evidence for a branch cut in the zeta function imply- 
ing slower-than-exponential decay of the trace. Finally we have identified a dominant 
timescale which is much longer than the shortest periodic orbit of the system. 



5 Concluding remarks 



In this paper we have discussed relations between periodic orbits, the resonance spectrum 
and the dynamical behaviour of a system. We have suggested a close relationship between 
the classical zeta function Z{k) and the corresponding function Z{k) as obtained from the 
BER approximation. We have not pursued this relation very far for the hyperbola billiard. 
This system is peculiar having infinite phase space volume so that the derivation in section 
2.2 becomes dubious. To study these ideas in detail will require much more theoretical 
work as well as numerical studies of a range of systems. Ergodic billiards calling to be 
investigated in this respect are e.g. the Sinai billiard and the closed three disk billiard. 
For the Sinai billiard it is straightforward to show that the function p(A) should decay 
as suggesting that the trace behaves asymptotically behaves as tr{t) — > 1 — const /t. 
Note that for this systems eq. (^) is not valid; the marginally stable orbit causes problems. 

For the case of the hyperbola billiard we observed that the relevant times scale was 
much larger than the period of the shortest cycle. The dominant timescale for the Sinai 
billiard should be the mean free path between the disks bounces which is also much longer 
than the shortest cycle, for small disk radii. 

So far we have only been discussing billiards. The resctriction to billiards is not only a 
convenient one (periodic orbits are much more easily found in billiards). Smooth systems 



seem to differ from billiards even qualitatively 19, |2^. The reason is that ergodicity 
seems to be lost at the momemt the billiard walls are made soft. To model this, it is 
convenient to embed the billiard into a family of smooth Hamiltonians, e.g. the hyperbola 
billiard is obtained as the a ^ limit of the family of potentials V{x,y) = ^(x'^y^)^/". 
When a is increased (the walls get softer) the finite sequences terminate earlier and earlier, 
indeed the sequences get considerably shorter when a = 1 as compared to the billiard 
case a = [19|. Orbits are pruned off familywise (c.f. the definition of family in section 
4) through bifurcations. In ref. |15| it is demonstrated how the bifurcation of a family 
belonging to a simple sequence also involve composed orbits. In each such bifuraction 
cascade, some orbits get stable and ergodicity is lost. It does not seem sensible that 
there can be any parameter value a where no such stabilization occurs. The occurence 
of stable islands seems to be a general feature of smooth potentials. 

A very interesting application of the trace concerns correlations in quantum spectra. 
The reason for this coupling is the close resemblance between the classical zeta function 
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(^) and the semiclassical one |21] 




The /Xp's are the Maslov indices and Sp are the action integrals which for homogenous 
potentials, such as billiards, are proportional to the periods Tp times some power of 
the energy E. The zeros of this object are (in the semiclassical hmit) the quantum 
eigenvalues. The spectral density, as given by the logarithmic derivative of Zgc, yields 



the Gutzwiller trace formula |22] 



One such interesting correlation measure is the spectral rigidity As{L). It is defined 
as follows. Given the quantum spectra {Ei}, consider the spectral staircase function 
N(E) = J2i^i^ ~ ^i)- Now, make the best fit, by means of the least square method, 
by a linear function over L mean level spacings. Calculate the square deviation between 
N{E) and this linear function, average this quantity over an energy interval AE which 
contains many levels but is semiclassically small (AE is smaller than the energy scale 
given by the shortest periodic obit in the system), and you end up with the function 
A3(L). Under the assumption that the trace, tr{t), is exactly one when t > Tmin Berry 



-'max 



[23 1 showed that A3(L) is consistent with GOE for a range < L < L^ax where L^, 
corresponds to Tmin-, if the system is chaotic and time reversible. But for bound chaotic 
systems we have realized that the trace may be significantly different from one up to a 
timescale much bigger that Tmin- This means that there is a region in L < L^ax where 
A3(L) cannot exhibit universal behaviour and where we can perform semiclassical calcu- 
lations. This is really what quantum chaology is about, to relate the quantum behaviour 
to the classical by semiclassical methods. Indeed, such departure from universality for 



L < Lfnax lias been observed in a number of numerical calculations |24, 25, 26 1. More- 
over, if there is a contribution from a branch cut of the zeta function we expect the 
trace tr{t) not to approach unity exponentially but rather as a powerlaw. This extends 
further the dynamical interesting region in L. Indeed one can show that if the trace goes 
asymptotically as tr(t) — > 1 — c/t the universal (GOE) result can only be achieved in the 
deep asymptotic limit E ^ oo. 



Time scales in connection with the spectral rigidity has been discussed in ref. |25|. 

We have seen that the trace for the hyperbola billiard is closer to 0.5 than 1.0, for a 
considerable range in t . This means that the A3(L) should be closer to the GUE result 
than to the GOE (which would be the appropriate since the system is time reversible); 
this is exactly what is observed in a quantum mechanical calculation [^6| . 

A lot of the recent work in the area of quantum chaology has been occupied with the 
problem of relating individual eigenstates to the periodic orbits. Due to the exponential 
proliferation of cycles this can not be pursued very high up in the spectrum. This 
problem is in itself very interesting. But, if we want to use semiclassical methods to 
investigate the transition from classical to quantum mechanics, that is, if we are interested 
in (asymptotically) high states, then we are bound to understand the asymptotics of the 
set of periodic orbits. The problem of exploring the asymptotics of this set is well 



examplified in ref. |27|. It is the author's hope that this paper provides some new ideas 
in this connection. 
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Appendix 



We now demonstrate how to obtain a reccurence relation for the coefficients of the 
pansion of 

k=0 n=0 

Consider first the simpler problem of expanding the Euler product 

oo ,00 
fe=0 n=0 

We extract the first factor 

00 , 00 , 00 

n (1 + -ri) = (1 + i) n (1 + Tl) = (1 + E «n(A)(*/A)" 

fc=0 k=l n=0 

and obtain the recurrence formula 



which can be exactly solved 



A 2 



" nLi(i-A-^) ■ 

We now use the some technique for eq. (^) 

k=0 k=0 k=l \n=l / \n=l 

leading to the reccurence relation 



bn = E an^lblA 



n-1 

-I 



.... 6n = 1 . 

1=0 
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Figure captions 



Figure 1: The spectrum (a) of model A (o) and B (+) and its asymptotic expression (□) 
and (b) the trace for model A (full line) and B (dashed) . The asymptotic expression has 
been used to calculate its discrete contribution ^ zeros to the trace. 

Figure 2: When Z(k) has a branch cut along the positive imaginary axis the Fourier 
transform of the logaritmic derivative has to be integrated along ^branch- 

Figure 3: The trace for model C 

Figure 4: p(A) (a), the spectrum (b), and the trace (c) for model D with different per- 
turbation strengths e. 

Figure 5: Members of some simple sequences in the hyperbola billiard 

Figure 6: Members of some composed sequences in the hyperbola billiard 

Figure 7: Graphical representation of the allowed cycles in some composed sequences 

Figure 8: Spectrum of the hyperbola billiard for various numbers of pseudo orbits N^ut- 

Figure 9: Location of the zero as a function of Ncut- 

Figure 10: The trace for the hyperbola billiard using direct summation (cq. 2) (long 
dashed), and as a sum over eigenvalues, with N^ut = 200 (short dashed), and Ncut = 460 
(fun hne). 

Figure 11: p(A) in the hyperbola billiard obtained by numerical simulation 
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Tables 





Jmax 




^ {jmax ) 


21(P 


3 


7.13= Ta 


19.65 


21210-'' 


16 


18.45= n 


527.9 


220^ 


51 


27.61= Tc 


922.9 


22W 


~640 


-92= Td 


~34000 


20^ 
IW 


oo 
oo 







Table 1: Some simple sequences of cycles. Invariants are given for last cycle in each 
sequence 
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